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Abstract 

In this paper we discuss the MHD flow of a second grade fluid, in particular we prove 
the existence and uniqueness of a weak solution of a time-dependent grade two fluid model 
in a two-dimensional Lipschitz domain. We follow the methodology of [3], i.e , we use 
a constructive method which can be adapted to the numerical analysis of finite-element 
schemes for solving this problem numerically. 
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1 Introduction 


A fluid of grade two is a non-Newtonian fluid of differential type introduced by Rivlin and 
Ericksen in [8]. An analysis in [1] shows that the equation of a fluid of grade two is given by 


2 . 

¥t 


[u — aAu) — vAu -|- {u — aAu)- S/uj — u V {u — aAu) = —Vp -|- / 
j 

div u = 0 


*Grupo de Matematica Aplicada, Dpto. de Ciencias Basicas, Facultad de Ciencias, Universidad del Bio-Bio, 
Campus Fernando May, Casilla 447, Chilian, Chile. E-mail: igor. kondrashuk@ubiobio. cl. I. K. was supported 
by Fondecyt (Chile) Grants Nos. 1040368, 1050512 and 1121030, by DIUBB (Chile) Grants Nos. 102609 and 
121909 GI/C-UBB. 

^Dpto de Matematicas, Facultad de Ciencias, Universidad de La Serena, La Serena, Chile. E-mail: 
enotteSuserena. cl. This author’s work was partially supported by project DIULS PR14151. 

^Dpto. de Matematica, Facultad de Ingenien'a, Universidad de Atacama, Copiapo, Chile. E-mail: 
mpobleteSmat .uda. cl. This author’s work was partially supported by Universidad de Atacama, project DIUDA- 
22256. 

^Grupo de Matematica Aplicada, Dpto. de Ciencias Basicas, Facultad de Ciencias, Universidad del Bio-Bi'o, 
Campus Fernando May, Casilla 447, Chilian, Chile. E-mail: markoSueubiobio. cl. This work was partially 
supported by project MTM2012-32325, Spain, Grant 1120260, Fondecyt-Chile and 121909 GI/C-UBB. 


1 



where a > 0 is a constant of material, z/ > 0 is the viscosity of the fluid, u is the velocity field, 
and p is pressure. For a = 0 the classical Navier-Stokes equation is obtained. 

On the other hand, in several situations the motion of incompressible electrical conducting 
fluid can be modeled by the magnetohydrodynamic equation, which correspond to the Navier- 
Stokes equations coupled with the Maxwell equations. In presence of a free motion of heavy 
ions, not directly due to the electrical field (see Schluter [4] and Pikelner [3]), the MHD equation 
can be reduced to 

^- Au + uVu——hVh = f - V{p* + —h^) 

Pm Pm Pm, 2 

dh 1 

-- Ah + u Vh — h ■ Vu = —grad oj (1) 

ot p,a 

div u = div h = 0 

with 

u\dn = h = 0 . ( 2 ) 

Here, u and h are respectively the unknown velocity and magnetic field; p* is the unknown 
hydrostatic pressure; oj is an unknown function related to the heavy ions (in such way that the 
density of electric current, jo, generated by this motion satisfies the relation rot jo = —aS/u) 
is the density of mass of the fluid (assumed to be a positive constant); p > 0 is the constant 
magnetic permeability of the medium; ct > 0 is the constant electric conductivity; z/ > 0 is the 
constant viscosity of the fluid; / is a given external force field. 

In the case the MHD equation coupled with the equation of an incompressible second grade 
fluid, the model can be write as 

d 

— {u — aAu) — uAu + curl {u — aAu) x u — (h ■ V) h = f — {p* -)- ft.^) 

C/ V 

o r 

— — Ah + {u ■ V) h — {h ■ V) u = —grad oj (^) 

div u = div h = 0 

with 

w Ian = h Ian = 0. (4) 

Note that when a = 0 we recover the model (1). 

One of the first mathematical results for this model type appears in [2], they prove the 
existence and uniqueness of solutions for a small time and global existence of solutions for small 
initial data in a conducting domain of based on the iterative scheme where discretization 
is performed in the spatial variables. In this paper we discuss the MHD flow of a second grade 
fluid, in particular we prove the existence and uniqueness of a weak solution of a time-dependent 
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grade two fluid model in a two-dimensional Lipschitz domain, where we follow the methodology 
of [3], i.e , we use semi-discretization in time and the work is in a domain of 


2 Preliminary results 

2.1 Notation 

Let (A:i,/c 2 ) denote a pair of non-negative integers, set \k\ = ki + k 2 and define the partial 
derivative by 

d V = —■. - j—. 

Then, for any non-negative integer m and number r > 1, recall the classical Sobolev space 

= {u G r(n);d^v G L'^{Q)y\k\ < m} , 

equipped with the seminorm 


and norm (for which it is a Banach space) 

1111 tv™.’’(O) = 




r 

E / 

d^v 

dx 





-1 1/r 


1 1/r 


E , 

0<|fc|<m 


n 




(n) 


with the usual extension when r = oo. When r = 2, this space is the Hilbert space The 

definitions of these spaces are extended straightforwardly to vectors, with the same notation, 
but with the following modification for the norms in the non-Hilbert case. Let u = (ui,U 2 ); 
then we set 


\u\ 


L’’(0) 


Un 


\u{x) 


1/r 


where ll'H denotes the Euclidean vector norm. 

For functions that vanish on the boundary, we define for any r > 1, 

Wo^’"(L!) = {uGWi’^(L!);u|ao = 0} 

and recall Poincare’s inequality, there exists a constant V such that 

Vu G iLo(H), lbllL’’(0) — blHi(o) • 


(5) 


More generally, recall the inequalities of Sobolev embeddings in two dimension, for each r G 
[2, oo) , there exists a constant Sr such that 


Vu G iLg(H), II'^IIl’’( 0) ^ ‘5r 


( 6 ) 
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The case r = oo is excluded and is replaced by, for any r > 2 there exists a constant such 
that 

VnGlTQ’(n), < Adr • C^) 

Owing to (5), we use the seminorm as a norm on Hq{Q) and we use it to define the 

norm of the dual space of Hq{Q,) : 




sup 

veH^in) 


{f,v) 


In addition to the norm, it will be convenient to define the following norm with the parameter 
a : 


V 



In the following, we denote by || • || the L?' norm. 

We shall also use the standard space for incompressible flow; 


H (div; 0) = {u G div v G 


H (curl; 17) = {u G curlu G L^(0)} 


y = {u G 7^0(0)^; divu = OinO} 

V^ = {ve Hl{nf-,yw G V,{Vv,Vw) = 0} 
Lq{^}) = {u G T^(n); qdx = O} 

and the space transport: 

X, = {fGL^ (0);u-V/gL2 (17)}, 
where u is a given velocity in . 


2.2 Auxiliary theoretical results 

To analize, we shall use the following results. The first theorem concerns the divergence operator 
in any dimension d. Its proof can be found for instance in Girault and Raviart [4]. 

Theorem 1 Let Q be a bounded Lipschitz-continuous domain of The divergence operator 
is an isomorphism from V-^ onto Lg (17) and there exists a constant (3 > 0 such that for all 
/ G Tq (17), there exists a unique v G V-^ satisfying 

dwv = f in 17 and ||^^|liyl(Q) < ^ ||/|| • 

The second result concerns the regularity of the Stokes operator in two dimensions, see [6]. 
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Theorem 2 Let Ll be a bounded polygon in the plane. 

1. For each r g]1,4/3[, the Stokes operator is an isomorphism from 

X onto 

i.e. for each f G L^{Q)‘^, there exists a constant Cr and a unique pair 
{u,p) G X n LliLi)] 

such that 

—vAu + Vp = f, divri = 0 in ri, ri = 0on dCl, 

and 

I^lvy2,r(0) + ll/llLr-l-Q) . 

2. If in addition, Ft is a convex polygon, then the Stokes operator is an isomorphism from 

nV X [i7^(r2) n io(^)] Furthermore, there exists a real number 

r > 2, depending on the largest inner angle of dll such that for all t G [2,r], the Stokes 
operator is an isomorphism from {W‘^’^{FL)Y FV x fl Lg(0)] onto L*{Fl)'^. 

The next result concerns the unique solvability of the steady transport equation in any 
dimension d, see [5]. 

Theorem 3 Let F be a bounded Lipschitz-continuous domain ofW^ and let u be a given velocity 
in V. 

1. For every f in L?‘{F) and every constant 7 > 0, the transport equation 

z + 7U • Vz = f in 12, 
has a unique solution z G and 

Ikll < ll/ll (8) 

2. The following Green’s formula holds: 

yz,9eXu, {u-vz,e) = -{u-xe,z). (9) 

Finally, the last result establishes compact embeddings in space and time. Its, proof, due to 
Simon, see [9]. 

Theorem 4 (Simon) Let X, E, Y be three Banach spaces with continuous embeddings: X C 
E CY, the imbedding of X into E being compact. Then for any number q G [l,oo], the space 

|uGL''(0,r;X);|^ GLi(0,r;y)| (10) 

is compactly imbedded in L^{0,T] E). 
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2.3 Formulation of the problem 


Let [0, T] be a time interval for some positive time T, let be an domain in two dimensions, with 
a Lipschitz-continuous boundary and let n denote the unit normal to d^l, pointing outside 
n. Let f G L^(0, T; iL(curl; 17)), the initial velocities UQ,hQ G V with curl (mq — oArio) £ 
and we expect the velocity u G L°°{0,T-,V) with du/dt G L‘^{0,T;V), the magnetic field 
h G L°°{0, t; V) with dh/dt G L^(0, T ; V), and the pressures p,uj ^ L^(0, T ; Lg(17)). 

The system (3) can be rewritten by introducing the auxiliary variable z = curl [u — aAu), 
as 

d 

— {u — aAu) — uAu + z X u — {h ■ V)h = / — V [p* + h^) 

3h 

—- Ah + {u ■ V) h — {h ■ V) u = —gradcj 


div u = div h = 0 


( 11 ) 


u = h = 0 on X (0, T) 


ri(0) = uq; h{0) = ho in 17 

Taking the acurl in the first of the above equations, it can be written as: 

a— + uz + a{u -V) z = z/curlri + acurl (h ■V)h + acurl/ — acurlV (p* + h^) (12) 

ot 

where we have used the fact that curl (z x u) = u ■'Vz, valid in two dimensions. Considering 
the above equation, we can rewrite the system (11) as follows: 
d 


— (w — aAu) — vAu + z X u = f + {h ■ V) h — V [p* + h?) 
— Ah + {u • V) h — {h ■ V) u = —gradcj 


dt 

dh 

'dt 


Oz 

a— + vz + a {u • S/) z — vcurlu = acurl (h • V) h + acurl/ 


(13) 


—acurlV [p* + h^) 


div u = div h = 0. 
Semi-discretization in time 

Let > 1 be an integer, define the time step k by 


k = 


T 

N 


and the subdivision points 7” = nk. For each n > 1, we approximate f (t^) by the average 
defined almost everywhere in 17 by 

/”(x) = ^ / f{x,s)ds. 

k Jtn-l 
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We set 


vP = uq, hP = Hq and = curl (uq — aAuo) . 

Then, our semi-discrete problem reads: Find sequences , (^")n>m(^")n>l,(P")n>l and 

such that u^,h'^ € V, € L^(fi), and G Lq{Q), solution of: 

I {u^+^ - u'^) - a\A - u^) - vAu^+^ + z^ x 

rC K 

= - V {p*^+^ + , 

I (h^+^ - h^) - Ah^+^ + ■ Vh^+^ - ■ Vw”+^ = -grada;"+\ (14) 

rC 

^ (z^+^ - z”) + iyz^+^ + au^+^ ■ = iycurlu^+^ + acurl/’^+i 

k 

Tacurl (h^+^ ■ Vh^+^) - acurlV (^p*^+^ + (h^+^f'j 

Now, we will make some estimates for u^,h^,z^, p* and w*. Multiplying the first Eq. of (14) by 
2ku^^^, the second Eq. of (14) by 2kh'^^^ and the third Eq. of (14) by 2kz'^^^ and observing 
that curlVE = 0, for any vector field E, we obtain 


2 — w*, — 2a (A (ri*+^ — w*) , — 2kv (Aw*+^, 

= 2k (/*+\ ri*+i) + 2k (/i*+^ • w*+^) , 

2 (/r*+i - h\ h^+^) - 2k {Ah^+^,h^+^) = 2k (/i*+^ • Vw*+\ , 

2 (E+i - zS ^ (E+i, E+i) = 2fe (curir+\ E+^) 

(curltr'+\ E+^) + 2fe (curl {h}+^ • V/r'+^) , z^+^) , 
where we used that fact that (ri*+^ • = 0 . 



7 



Proof; Multiplying the first Eq. of (14) by and the second Eq. of (14) by 2/1*“'“^, we 

obtain 

^ - a^A - 2v (A'i/*+\ 

= 2 (/*+\ ri*+i) + 2 , 

^ {h^+^ - h\ - 2 - 2 (^*+1 • Vri*+\ = 0 

where we should note that 

(z* X w*+^w*+l) = 0, (^V = 0, 

(grad a;*+i, = 0, {u^+^ ■ V/i*+^, = 0. 

Using the formula 

2(a - b,a) = ||a||^ - \\b\\^ + ||a - 6||^ , (17) 

that is true in any Hilbert space, and adding the above equations, adding from i = 0 to n — 1 
and making use the telescopic property, we have 


- ||ri”|P + - 
A:" "“^A: 


n—1 


2 1 
k 


n—1 


+ — I ri"'" — ri" 


n—1 


,*+i 




i=0 
n—1 


i=0 


+2^J^i|Vu«|| +2^||Vh'+'||' 


2=0 

n—1 


2=0 


< 


2CEll/ 


*+i|| IUd+1 


u 


+ T ll'^olla + T 


2=0 


k 


k 


Now taking into account that 


2C||/ 


*+i|| lk*+MI < 


(18) 


4^2^||y.m||2 1 ||^2+in2 

2 " 26 " 


<4c^|l|r«l|Uk||v„«f 


then from equation (18) we can write 

n—1 n—1 

2i/^ ||V-a*+^|f + 2 ^ ||V/j.*+^||^ < 


2=0 


2=0 


n—1 


2=0 




,2+11 


2(5 


1 


1 


+ tKII: + 7II^oII 


k 


k 


then, putting 6 = CI2v, we obtain 


n—1 


n—1 


■EI|V«-+‘|| + 25 ;||Vh'«|| 


2 II J,||2 , 1 ||_ ||2 , Mil, l|2 


2=0 


2=0 


vk 


ll/llL2(n,x]0,t’2[) + T ll'^olU + T ll^oll • 



On the other hand, to obtain estimates of the || , we multiply the second equation 

in (14) by 2Ah'^^^ ^ then we obtain (after applying the projection operator P) 

^ - Vh\Vh^+^) + ^ = -2 A/r*+^)+2 , 

then bounded each of terms, we have 


1 


k 


(V/i*+i - Vh\ V^*+^) = - ||V^*+^|| - - \\VhY + 1 ||V^*+^ - 


k 


k 


\ 2 {u^+^ ■Vh^+\Ah^+^)\ < 2 || ii *+^||^6 || Vh *+^||^3 || A / i *+^|| 

< 2 ||Vii*+^||^ \\Ah^+^\f^^ 

< 2Ce, ||Vw*+^||^ + 2^1 11^^' 


i+l\ 


\2{h^+^ ■Vu^+\Ah^+^)\ < 2||/i*+i||^^ ||Vw*+^|| ||^/i*+^|| 


I V2 II 4^*+1||3/2 

< 2CC„ ||Vw*+^ir l|V/i*+^||^ + 2C£2 > 


< 2C\\Vh^+^\f^\\Ah^+^\f 

,i+l||4 


|Vw*+^| 

|2 . „ „ II , , i-^-^ l|2 


where we use the estimate of interpolation ||^*'*'^| \l^<C . From above 

estimates and taking into account that ||Vti*^^|| is bounded, we have 

||Vfi*+^||^ + ||V/i*+^ - Vh^f + || ||2l/i*+^|f 

< A:(2C'^, +2CC£2)C||V/i*+^|f + A:(2ei + 2Ce2) 

then, there is k and e such that (for a sufficiently large N) 1 — k (2C'ej + 2CC£^) (7 = 1/2 and 
2k/jia — k (2ei + 2 Ce 2 ) = k/ixa, thus, from the above inequality we can write 

2 ||V^*+M|^ + ||V^*+^ - < ||V^*||^ , 

/ifj 

from which we get (using the lemma 3.14 pg. 131 in [11] with r] = Ji = ^ = 0) 

12 , ^ iIt-^. >-^,iii2 , k 11., ||2 


2 ||v/i*+^|f + ^ M|V^^+1 - vwf + ) < ||V/iof • 

i=i ^ 


□ 


Proposition 6 The sequence {vA)n>i and {hA)n>i satisfy the following uniform a priori esti¬ 
mates, for 1 <n < N 


\u'' 


n—1 

+ Ei 

n —1 




vd\ ^ ll/llL2(Qx(o,t")) + ll«o||^ + ||^of 




+ X] “ ^11 ^ ^ ll/llL 2 (Ox( 0 ,t")) + ll'“olla + ll^ol 


(19) 

( 20 ) 


i=0 
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Proof: Estimates (19) and (20) are derived by adding the first and second equations of (15) 
and using the formula (17), 

- \\h^\f + - h^f - 2a (A {u^+^ - u^) ,u^+^ - w*) 

-2a (A («*+i - u^),u^) - 2kv + 2k ||Vh*+^||^ = 2k (/*+\-a*+i) 

then adding from i = 0tof = n — 1 and making use of the telescopic property, and again using 
the formulae (17), we have 


W' 


2 

a 


n—1 


'^oWa + X] 

i=0 



2 


n—1 


n—1 


l^of + j;||/i'+'-/i'|r + 2aj;||V 


- u^) 


i=0 


i=0 


( 21 ) 


n—1 


+2t^l|vh'+'|p< 

2=0 2=0 


2+11 


n—1 


n—1 


now. 


drooping 2a ||V(++i -+)|| and 2A:X] ||Vh*+i|| we have 


i=0 


i=0 


U" 


n—1 

+ E 

2=0 


n—1 


I++1 -+| 




i=0 


- ^ ll/llL2(Q,x]0,t"[) + ll'“o||« + ll^oll^ 


n—1 


where ^ | ^ = IIZIIl2(o x]o t"[) • which we get the result. 


i=0 


□ 


Remark 7 From (21) we have 

n—1 n—1 ^2 

2o ^ !| V {i.‘+‘ - u‘) II" + 2t II V/."f < — + lltioC + ll'>olf ■ (22) 

2=0 2=0 

Estimate for z” 

Before obtaining an estimate for , we will show the following corollary 
Corollary 8 Give the sequence (h,"')^>^ we have the following uniform a priori estimates 
\\Ah^<C and ||V (^’^ • Vh’^)|| < C Vn G N. 
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Proof: Applying the projector P to the second equation of (14) we obtain 

= -P {u^+^ • + P {h^+^ • - ip - /i”) 


||A/i'^+i|| < • V/i^+i|| + \\h'^+^ • Vri^+i|| + i \\h'^+^ - h" 

thus, each term can be estimated as follows: 


•Vh’^+^ll < 


< C\\Vu^+^\\\\Vh^+^f^\\Ah^+^f^ 

< Csi + El forei > 0 small 


here was used L® and a result of interpolation. 


/jn+l . vtt^+l|| < ||/i«+l| 


— 11”' llL“(0) 


llL°°(0) 

n+l||l/2 II . 


< C,,+e2\\Ah^+^\\. 


Finally, substituting in (23) we obtain 


(1 - £1 - £ 2 ) \\Ah^^^\\ < Ce, +C,, + C 
and considering (1 — £i — £ 2 ) > 0 we obtain 

||A/i'^+i|| < C, 

where the constant C is generic. 

Now, taking into account the equation (24) and usual estimates, we have 

< Cl ||V/i”+^ • V/i^+^ll + C2 ■ V^h^+^\\ 

< C5 ||A/i'^+i|| ||Ah^+i|| AC'e/^ ||A/i'^+i|| 
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another consequence of (24) is ft."" • < C. 


□ 


Proposition 9 The sequence (-z"')n>i satisfy the following uniform a priori estimates 


n—1 

|2 , II i+1 




i=0 

(C^c ...2 y I, „2 1 II, II 

- ^ ll^olL + - ll^oll 

4 —— l|curl/||^2(f2x(o,t")) + • 

Proof: Multiplying the third equation of (14) by we get 


2( 


z ' — z , z 


I s 2/cz^ / 

^) +-(r' 

^ a ^ 


1) = ^(curln*+i,z*+i) 


(25) 


+2A: (curl + 2k (curl ■ Vh^+^) , z^+^) 

then, using the formula (17), we obtain 


7*+l I 


2ki' 




2kn 




*+i| 


<-||Vn*+i|| ||z*+i|| + 2k ||curl/*+i|| ||z 


a " 

*+i|| IU*+i| 


a 


+2A:||V(/i*+^ • V/i*+i)|| ||z*+i 
2ki' 


< 


a 


(Cei ||Vn*+^f+ e? +2A:C£2 ||curlf+i 


+2A:e2 + 2kC\\z 


i+l| 


2kv 11^ , 
<-I Vtt 


m|| 2 ^:^|U*+i|| 2 ^ 2 —||curir+H 

2aei " a 2e2 


+te2 ||^*^^ir + — + fcealiz 


\2^k 

£3 

Now adding to i = 0 to n — 1 we obtain 

n—1 


. 2+11 


E 

i=0 


j+i| 


- u; 


n—1 

+E 

i=0 


U*+i -z*| 


(2kv ku 7 ^ \ 

+- £1 - k £2 - kCss > 

V a a / 

X ■' i=0 

si|Eiiv«-ir+"fi|c..o 

2=0 2 = 0 


d+i| 


*+l||2 + C-. 


£3 
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n—1 


k -v ^ k k T T 

Where the term — y C = —Cn < — C— = C —. Now considering ei = 1 and £2 = £3 = ^^/2c 

£3 ^ £3 £3 A: £3 

2=(J 

the above equation can be written as 


n —1 

E 

i=0 


.*+ll 


n—1 


- U: 


n —1 

+ Ei 

i=0 
n—1 


E+i - E 


S ^E l|vt.‘«ir + ^E 

i=0 2=0 


consequently, 


n—1 

+ EI 


- Ell < 


i=0 


ll/llL 2 (O,x] 0 ,t"[) + ^ ll^olla + - ll^ol 


□ 


Proposition 10 Let 


Cz = sup IIE 
0<n<Ar-l 


The sequences — u^)/and — h^)/, satisfy the following uniform a 

priori estimates: 


1 


n—1 

'^2k 

i=0 


|n*+i - w* 


2 ^ (Cii/ + cvl) 


2 . 


2a 


f ^2 lllllL2(o,x]0,t"[) 


+ ~ ll'^^O 

V 


2 

a 



+ cl ll/lli2(Qx]o,t"[) + C>T, 


n—1 - 

Es 


i=0 




2 ^ {Ciu + C,Slf (C^C 


2a 




+“ II'^^oIIq + - ll^oll^^ + Cl ||/|li2(f2,x]o,t"[) + C>T. 

Proof: Multiplying the first Eq. (14) by — n*) and the second Eq. 

(/i*+i - h^) , we obtain 

\ - vf) - (n*+i - u\ - vf) - v (Aw*+E - vf) 

+ (z* X - u*) = - u*) + (^*+1 • - u*) , 

i - h^) - (A/i*+E - hf) + - /i*) 

- - hf) = 0, 


of (14) by 


(26) 
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then adding the above equations and using that 


Z* X — uM < CzS'f u 


'2 IL,i+l| 


7#*+^ - 71* 


we have 


- ||w*+i -7x*||^ + i ||/i*+i -/i*| 


1 

k 


“ ' A: 

2 IU,i+i||2 


<(Ciz. + a5|) — im 


a 


,i+l 




+ 


+ 


*^ 2^2 II ^ i + l ||2 , J ^||^_ i +1 ^_ i ||2 ^ g|fcg3 ||^i+l y^i+l||2 


\\f 


2ke3 


u 


i+l 


,, H- \\u'’^-u''\\ + 

II ^ 2e2 " II 2 


+-^ ||/i*+i - + C^k— ||-u*+i • Vh^+^\f 

2A^£‘4 2 




|/i*+i • Vw 


i+l I 


+ 


2 keQ 


\h^+^-h^\f 


then put El = k, 62 = 2k, 83 = 2 ,64 = 2, £5 = 4 and £7 = 4, of the above equation can be written 


— ||w*+i - - Ai* 

2A: II 11“ 2k II 


2t2 ^ IL.i+l|'2 


<(Cii/ + a5|) — Iltr 


ii(o) + Cik ||r+i|| + Csk ||/i'+i • (27) 


+C|A: 11 ^/ 1 *+^ 11 + 2C'|A: ||7i*+i • Vh^+^\\ + 2Clk ||/i*+i • V7x*+i||^ 
On the other hand, recalling the Corollary 7, we have 

|2 , , II . , i-Ll I|2 


|/i*+i-V/i*+^|| <di, ||7l/i*+i|| <d2 


\U' 


i+l 


VAi*+i||^ < ds, 11^*+^ • Vu*+i||^ < di 


and summing from i = 0 to n — 1 in (27), we obtain 


n—1 - n—1 - 


i=0 
n—1 


i=0 


<E(C: 

^=0 

n—1 

+ Y,Dk, 


2i2 


n—1 


2 a 




i+l I 




+c'IEHI/* 


i+l I 


i=0 


i=0 
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Proposition 11 The sequence (p"')„>i and {u}^)n>i satisfy the following uniform a priori esti¬ 
mates: 

n—1 

i=0 


'{Ciu + C.SjfL 
a 




+ [4C|C + 2C|L] ll/f + [D + 2LD]T, l<n<N. 

Proof; Let vi,V 2 G V-^ = {v G Hq (11) ;Vtr G V, (Vv, Vtr) = O} such that divui = and 
gradu2 = Then by multiplying the first and second eq.(14) by vi and V 2 respectively, and 
adding the results, we have 

i {u^+^-u\vi) + (z* X u^+\vi) + 

- (^*+1 • Vh^+\vi) + j - h\v2) + • Vh^+\v2) 

K 

- ■ Vn*+i,U2) + (a;*+i,u;*+i) = (/*+\ui) . 
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where, we consider that vi,V2 G V'^. Thus, we can write 

k + k < Ciy ^ lluif 

^ lluill^i + kcfj- ||/r*+i • V/r*+i|f 

+A:C|| ||ri*+i • ^ ||/i*+i • Vw*+i 

Now, we considering that M- L^, i.e., 

||t^i||<C||Vr;i||=C||p*+i|| 


and 

||t'2|| ^CllVuall =C||u;'+^||. 

From which, we obtain 

k ||p*+^|p + k ||w*+^|f < Cl— ||tr*+i - 

+tc?s|| ||t.‘+i^. + T ||p-+'||" + )cc|| ||w+‘. v/.‘+'|p 

+g ||p'«f + c|| II v+‘ - h'lp + TIP 

+tc|^ ||u‘+' . Vh‘+'|p + ^ ||w‘+"|p + kC'i^ ||W+‘ ■ Vu‘+"|p 

2 2u^ 2 


Ck 


,*+ii 


’'^7 II j.i+ii|2 , CA: II j_,_i||2 


+kcri\\f^r+^¥ 


2 " ' 267 

then, taking ( 5 i = AC/k, 62 = A, 5 ^ = 67 = AC and ^4 = 2 C/A, ^5 = 
to n — 1 and multiplying by 2, we have 

n—1 n—1 


A:||p*+^||^ + A||rc*+^ 


i=0 


1 


S‘‘TcEr 


i=0 


u — u 


AC, addying from i = 0 


+4C25| k\\u^Y\l^ + 2ClcYv - hf + ACiC \f\h^n,]o,Tl) 

i=0 i=0 ^ 

n—1 

+ Y {^ClCdi + 4C|Cd3 + 4C|Cd4) k, 

i=0 
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thus, put D = AC^Cdi + AC^Cds + AC^Cd^, noting that n < T/k and using the inequalities 
(28) and (16) , we have 


n—1 


k\\p^+^\f + k\\oj^+^\f 


< 


i=0 


< 


{Cy^ + C.SifL , 


+ iC^Si 


a 


( C^C 1 1 

^l/li^(Ox]0,t»[) + -KI|^ + -||^0f 


+ + 2C'|ii] \f\L^(nx]o,T[) + [D + 2LD] T. 

where L = max {dCfC, dCfC} 


□ 


2.4 Existence of solutions 

Here, it is convenient to transform the sequence {u'^), {h"^), (p^), (<^”) and (z^) into functions. 
Since (h^) and (z”) need to be differentiated, we define the piecewise linear functions in 
time; 

■j- _ 

Vte[r,r+1], Uk{t)=u^ + — —(m^+I-w^), 0<n<N-l 

k 

■f _ -f-^ 

VtG[r,r+i], hk{t) = + ^—- h'^) , 0<n<N-l 

k 

-j- _ 

yt G [r, Zkit) = z^ + —— -z^), 0<n<N-l. 

K 

Next, in view of the other terms in (14), we define the step functions: 
yt G [t"^, and 0 < n < — 1; 

/fc(t) = r+i, Wk{t) = u-+\ g^{t) = h^+\ p^{t)=p^+\ 

Uk{t) = Cfc(t) = Afc(t) = z'^. 

Then we have the following convergences. 

Proposition 12 The exist functions u,h £ L°°{0, T ; V) with du/dt, dh/dt G T^(0, T;V),p,u} G 
L^(0, T; Lq(11)) and z G L°°(0, T; L^(n)) such that a subsequence of k, still denoted by k, satis- 
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fies: 


lim uj. = lim = u weakly * in L°°{0, T;V), 
fc—>-0 fc—>0 


lim hk = lim g^ = h weakly * in L°°(0, T;V), 
fc—>-0 fc—>-0 


lim Zk = lim = lim = z weakly * in L°°(0,T; L‘^{Q,)), 
fe—>-0 fc—>0 fe—>-0 


limpi~=p weakly in L^(0,T; LQ(r2)), 
fe—>-0 


lim Wfc = w weakly in L^(0, T; LQ(ri)), 
fc—>-0 


, d d 
lim —Ufc = —u 
fe-i-0 at at 


lim —h 

k->-0 dt 


k = 



weakly in L^(0,T;F), 
weakly in L‘^{0,T;V). 


Furthermore, 

lim Wfc = lim Wfc = u strongly in L°°(0, T; L^(n)^), 
k^O k^O 


lim/ifc = limgr^ = h strongly in L°°(0, T; L^(n)^) 
fc-lO fc-iO 


(29) 


Proof: Due to the uniform estimates given in Propositions 5 -10, we can extract a subse¬ 
quence (still denoted by k) such that; 


lim tifc = ri; lim = h, weakly * in L°°(0, T; D), 
fc—lO fe—>-0 


\mi Zk = z weakly * in L°°(0, T; L^(D)), 
fc-iO 


limpfc=p; lim uifc = tc weakly in L^(0, T; Lq(D)), 
fc-iO fc-iO 


lim —Wfc = —w; \\m—hk =—h weakly in L2(0, T; D), 

fc-io at at fc-j-o at at 


lim rufc = w; lim Qk = 9 weakly * in L°°(0,T; V), 
fc-iO fc-lO 


lim lim Afc = A weakly in L^(0, T; L^(D)). 

fc-iO fc-lO 
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As far as the function C, and A are concerned, observe that 


1_ -j- 


VtG[r,r+l], g^{t) = hu{t) = 


k 

- t 
k 


\u 


— u^) , 0 < n < N — 1 

-h^), 0<n<N-l 


- t 


VtG[r,r+i], Ckit)-Zkit) = --- (z^+^-z^), 0<n<N-l 


- t 


\/t e Xk{t) - Zk{t) = - - -(2;"-+^ - z”) , 0<n<N-l. 


Therefore 


N-l 


W'Wk '^k\\L'2{0,T,V) ~ o 


n=0 


„,n+l _ n 2 

H^(n) 


N-l 


\\gk-hkf = -Y.\\h^^^-h^^ 


n=0 


\m{n) 


k 

IlCfe — ^A:lll,2(Qx]0,r[) = W^k — ^fc|lL2(Qx]0,r[) = T “ ^""11 ' 


n=0 


(30) 


Then, using the estimates (19),(20),(25) and the uniqueness of the limit, we have w = u,g = h 
and C = A = z. It remains to prove the strong convergence (29). In view of (30), it suffices to 
prove the strong convergence of Uk and hk- Note that, (tt^) and (hk) are bonded uniformly in 
the space 

G L\o,T-,H^infy, ^ G L\0,T-,L\nf)^ , 

and as the imbedding of into L^(II) is compact, the Simon’s theorem implies that Uk and 

hk converges strongly to u and h respectively in L^(0, T; L^(II)^). □ 


Theorem 13 Let LI he a bounded Lipschitz-continuous domain in two dimensions. Then for 

any a > 0,u > 0, / G L^{0,T; H{cuil;Ll)) and uq^Hq G V with curl (wq — oAwo) G L‘^{Ll), 

d'Ui dhy 

problem (11) has at least one solution u,h G L°°(0,T;y) with € L'^{0,T-,V) and 


p,cuGL2(0,r;L2 (II)). 
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Proof: Let /c be a subsequence satisfying the convergences of above Proposition. It is easy 
to check that the functions u^, Zk^Pk-,^k-, Wk, Qki Ck and satisfy the following formulations; 


VuGiLi(^^),V(/pGCO([0,r]), 


d 


^^Uk{t),vj + a ( —Vwfc(t), Vu 


d 


+v (Vwkit), Vv) + (Ak(t) X Wk(t),v) - {pk{t),diyv)] if{t)dt 


(31) 


+ / {9k{t) ■'^v,gk{t))p{t)dt = {fk{t),v)ip{t)dt, 

Jo Jo 


ygeH^{n),ycj,GC^{%T]), 


rT r 


^hk{t),g ) + (Vsifc(t), Vgr) 


-{wk{t) ■ Vg, g^it)) + (gkit) ■ Vg, Wk{t)) - {w^, divg)](p{t)dt = 0, 

ye G ipi’4(L!),vV’ G ci([o,r]) with V’(t) = o, 

-a j {zkit),9)-^i;{t)dt + j [v {Ck{t),0) - a {wk{t) ■ Ve,Ckit))]'il^{t)dt 
— [z^, e) V'(O) = f [v {curlwk{t),9) + a (curl/^(t), 0)] ^{t)dt 


(32) 


(33) 


+ / a{cml{gk{t) -Vgkit)) ,9)'ilj{t)dt, 
Jo 


where we note that 


= I - “") 

Note that, the weak convergences of the proposition above, imply the convergences of all the 
linear terms in (31),(32) and (33) and the terms involving f also converge, from standard 
integration results. Thus, it suffices to check the convergence of the non-linear terms. Then, for 
all indices i and j, I < i,j < 2, 


d 


d 


~ dt 


t-f 


u 




and 


lim {wk)iVj(p = UiVjip; 
k—>-0 


lim (gk)i gjct) 


higjfj) strongly in L^(nx]0,T[) 


lim Xk = z weakly in (11 x ]0,T[), 
k^O 
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then, we have 


Similary 


Therefore 


rT i-T 

lim / {Xk{t) X Wk{t),v) ip{t)dt = / {z{t) x u{t),v) (p{t)dt, 

Jo Jo 

lim / (fl'fcW • = / {h{t) ■Vv,h{t))<f{t)dt, 

7 o Jo 

lim / (■mfc(t) • Vg, g^it)) (p{t)dt = [ {u{t) ■ Vg, h{t)) (p{t)dt, 
Jo Jo 

lim / (gkit) ■X/g,Wkit))(j){t)dt= [ {h{t) ■ Vg,u{t)) (j){t)dt. 

Jo Jo 

lim {wk ■ V0) V’ = (w • V0) V’ strongly in L^(nx]0, T[), 
k^O 

lim {g;. ■ Vcurl0)'0 = {h ■ Vcurl0)'0 strongly in L^(nx]0,T[) 
k^O 

lim / {w}^ ■V6Xk)'4^{t)dt = [ {u ■ VO, z{t)) ^{t)dt, 

Jo Jo 


lim [ (gj^ ■ VcuTl9,gk) ^{t)dt = [ {h ■'Vcurl6,h)'ilj{t)dt. 

Jo Jo 

Hence we can pass to the limit in (31), (32) and (33) and we obtain 


^v€H^{n),Vip€C^{[0,T]), 


dt' 


0 

u{t),v ) + a { — Vw(t), Vv 


+v (Vri(t), Vv) + {z(t) X u(t),v) — (p(t), divi))] ip{t)dt 

+ [ {h{t)-Vv,h{t))ip{t)dt= [ {f{t),v)ip{t)dt, 

Jo Jo 


^g€H^in),^<jGC%[0,T]), 


A 

w 


Ht),9 + (V/j.(t), Vfir) 


- {u{t) ■ Vg, h{t)) + {h{t) ■ Vg, u{t)) - {uj, divgr)] 4>{t)dt = 0, 
V0g1H1'^(O),VV^gC1([O,T]) with v^(r) = o, 

—a I {z{t),9) ■^'ip{t)dt + I [v {z{t),9) — a{u{t) ■V9,z{t))]^lJ{t)dt 
Jo d)t Jq 

— 0) V'(O) = f [v (curlri(t), 0) + a (curl/(t), 0)] 'i(j{t)dt 


+ / a {curl {h{t) ■Vh{t)) ,9)'ilj{t)dt. 
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By choosing v,g £ 'D{Q)‘^, (p, cj) and € ^^(]0, r[) and 9 £ we easily recover (13). It remain 

to recover the initial data, for this note for any g £ L^(n)^ and any (/> £ satisfying 

t — 

4>{T) = 0, and used the formula hk{t) = H- - — — h"') we have 

K 

p'Y' / f) \ P) 

Jo = {hkit),g)—(p{t)dt-{h^,g)p{0) 

Passing to the limit in the equality above, we have 

/ p\ \ P) 

(Hi), 9) V?(0). 

On the other hand, we have 

F) r'^ Pi 

Jo {Kt),9)-Q^Pt)dt-{h{9i),g)p{f)) 

where we conclude that = h,(0), similarly we obtain vP = ii(0) and = z(0). □ 

With respect to the uniqueness it is possible to show an analogous to the [3]. In fact, we 
have 

Theorem 14 Assume that Q is a convex polygon. Then for any a > 0,u > 0, / in L^(0, T; i7(curl; II)) 
and uo £ V,h £ V with curl(rio “ aAuo) £ problem (3)-(4) has exaetly one solution 

{u, h,p,uj) £ W. 
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